AMENDMENTS TO THE CLAIMS 

1-18. (Cancelled) 



19. (New) A method for modelling, analyzing and/or synthesizing, a windowed 
signal, comprising computing simultaneously the frequencies and complex amplitudes 
from the signal using a nonlinear least squares method, whereby the computational 
complexity is reduced by taking into account the bandlimited property of the window 
resulting in band-diagonal system matrices for the computation of the amplitudes and 
frequency optimization step. 



20. (New) The method according to claim 19 using a stationary nonharmonic signal 
model according to (Eq. (2)): 

K-l 



w n A * exp(-2niuj k ^) 



k=0 



(2) 



which is a model with K stationary components where each component is characterized 
by its complex amplitude A k and frequency cu*, where w„ is the window of claim 1 ; or 



an harmonic signal model according to (Eq. (3)): 
5-1 s k -\ 

k=Q p=0 



(3) 



which is a model with S quasi-periodic stationary sound sources with a fundamental 
frequency co kl each consisting of S k sinusoidal components with frequencies that are 
integer multiples of co ki in which the complex amplitude of the pth component of the kth 
source is denoted A k>p , and where w„ is the window of claim 1. 



21. (New) The method according to claim 19 using a nonstationary nonharmonic 
model according to (Eq. (4)): 
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which is a model with K nonstationary sinusoidal components which have independent 
frequencies co k , in which the amplitudes A k>p denote the p4h order of the £-th sinusoid, 
and where w„ is the window of claim 1 . 



22. (New) The method according to claim 20, comprising the computation of the 
spectrum as a linear combination of the frequency responses of the window according to 
(Eq. (11)): 



X m = A kW(m + u>k) 



(11) 



k=0 



for the stationary nonharmonic model, 
or(Eq.(12)): 

S-l S k -i 

*™ = EE^( m +w) ( 12 ) 

fc=0 p=0 

for the harmonic model, 

where the fourier transform of a complex signal results in a spectrum Xm, where W(m) 
denotes the discrete time fourier transform of w n and whereby only the main lobes of the 
responses are computed by using look-up tables. 

23. (New) The method according to claim 21, comprising the computation of the 
spectrum as a linear combination of the frequency responses of the window according to 
(Eq.(13)): 
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(13) 



for the nonstationary model, where the fourier transform of a complex signal results in a 
spectrum 4, where W(m) denotes the discrete time fourier transform of w n whereby 
only the main lobes of the responses are computed by using look-up tables. 
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24. (New) The method according to claim 20, comprising the step of computing the 
stationary complex amplitudes, by solving the equations (Eq. (19)): 



(19) 



' B l,l 


B 1 ' 2 ' 




A r 




' C 1 ' 


B 21 


B 2 - 2 




A' 




c 2 



where 
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using (Eq. (20)): 

BH = i^(n^ + ^)) + ^(m-^)) 
fi/i = -^(n^+^))-^(n^-^)) 
^ = ~9(y(w*+wi)) + W(w*-wi)) 
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(20) 



m=0 



such that only the elements around the diagonal of B are taken into account, whereby a 
shifted form B j s computed containing only D diagonal bands of B according to (Eq. 
(27)): 
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- B l,l+k-D 

= B?? +k -D (27) 

and Eq. (20), whereby the computation of the Eq. (20) requires the computation of the 
frequency response of the window and the square window denoted by W(m) and Y(m) 

respectively, and solving equation given by Eq. (19) directly from B and C in (Eq. (28)): 

A r = SOLVE(B^~\C 1 ) 

A} = SOLVE(B?*,C 2 ) (28) 
by an adapted gaussian elimination procedure. 

25. (New) The method according to claim 20, further comprising the step of 
optimizing the frequencies for the stationary nonharmonic model by solving the equation 
(Eq. (34)): 

HAcj = h (34) 
using (Eq. (42)): 

h = -~»(i4, ^ -*»))) 

\ m=0 / 

H lk = X(A k A l Y'\u k + Cu l ))-U(A k A* l Y''(w k -u> l )) 
-Ajiu^S (A RnWfa - m) J + 5 fci A 2 

\ m=0 / 

(42) 

such that only elements around the diagonal of H are taken into account, whereby a 
shifted form H is computed containing only D diagonal bands according to (Eq. (36)): 

Hik = Hi,i+k-D (36) 

and Eq. (42), whereby the gradient h is computed from the residual spectrum R mi where 
= x m - x m denotes the spectrum of the residual r ny and from amplitude A/ and 
frequencies cc>u and requires the computation of derivative of the frequency response of 
the window W'(m) t whereby the first term of H requires the computation of the second 
derivative of the frequency response of the square window denoted Y"(m), whereby the 
second term of H is computed from the residual spectrum R mt amplitude Af and 
frequencies cou and requires the computation of the second derivative of the frequency 
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response W"(m) % whereby the parameter X\ allows to switch between different 
optimization methods and the parameter \ 2 regularizes the system matrix, and 

computing the optimization step by solving the system of equations directly on Hand h 
according to (Eq. (37)): 

MD = SOLVE(%h) (37) 
by an adapted gaussian elimination procedure. 

26. (New) The method according to claim 20, further comprising the step of 
optimizing the frequencies for the harmonic signal model, by computing the optimization 
step solving (Eq. (48)): 
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with 
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whereby the gradient h is computed from the residual spectrum R m% where 
Rm = x m - x m denotes the spectrum of the residual r n and W\m), and from amplitude A\ 
and frequencies co/, and requires the computation of derivative of the frequency 
response of the window W(m) M whereby the first term of H requires the computation of 
the second derivative of the frequency response of the square window denoted Y"(m) t 
whereby the second term of H is computed from the residual spectrum R mi amplitude A/ 
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and frequencies co/, and requires the computation of the second derivative of the 
frequency response W"(m) } whereby the parameter X\ allows to switch between different 
optimization methods and the parameter X 2 regularizes the system matrix. 



27. (New) The method according to claim 19, further comprising the step of 
computing the polynomial complex amplitudes by solving the equation (Eq. (55)): 



(55) 
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using (Eq. (63)): 
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(63) 



such that only the elements around the diagonal of B are taken into account, whereby a 

shifted form B is computed containing only PD diagonal bands of B according to (Eq. 
(64)): 



& lP+q,kP+p — £>/d , „ ;d , „ , i, d i ^ ... / n i i\d , i — & 
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and Eq. (63), whereby the computation is required of the frequency response of the 
square window and its derivatives d^ Y ^ whereby the computation is required of the 
frequency response of the window and its derivatives ~^ w ( m \ and solving the equation 
given by Eq. (55) directly from B and C by an adapted gaussian elimination procedure. 
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28. (New) The method according to claim 24 or 27, comprising a preprocessing step 
which comprises: 

sorting the frequencies to obtain a band diagonal matrix Z), determining the number of 
diagonal bands D being defined as the largest k - I for which ~P < ^k-^i < P t where 
u) k and cu/ denote two frequency values and j3 the width of the main lobe of the 
frequency response of the window. 

29. (New) The method according to claim 19, further comprising the step of 
computing instantaneous frequencies and instantaneous amplitudes according to (Eq. 

(69)):, 



*t(n) = 



'p-i \ 6 /p-i 



E • 4 ;,»<- 2 * =¥■>* + E A u-?'i^y 



p=0 



E; = o4, P (-27r^p) 



$ fc (n) = 27r^ fc ^p + i arctan . ,pV (69) 



whereby the instantaneous frequency can be used as a frequency estimate for the next 
iteration as expressed in (Eq. (73)): 

A r A* — A r 
(r+l) _ (r) /m Vm ^,0^,1 ( 7 o\ 

u k - ^ - Vn) — j~~2 — — \' 6 ) 



30. (New) The method according to claim 19, further comprising the step of 
computing damping factor according to (Eq. (78)): 

Ar at i Ai Ai 
, x ^k O^A: 1 ' ^k fl^fc 1 ,„^ N 

Pfc ~ - ^ -^7^5 : M 2 M ( 7 8) 

/U i Ar 

in case that the amplitudes are exponentially damped. 



31. (New) The method according to claim 19, where a scaled frequency response is 
used for the analysis of a zero padded window, where W M (m - m 0 ) denotes the 
frequency response of the window of length M and w m(™> ~ n o) = W A, (# m - m 0 ) the zero 
padded version of this window up to a length N t and ™M'( n ~ n o) = ^^ M (^n-m 0 ) t h e 
inverse transform of the truncated spectrum to a length N' reducing the window length to 
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M f — M ft/ 7 

1V1 n , resulting in a scaled and zero padded version of the window by computing 
the inverse transform of the scaled frequency response yielding (Eq. (1)): 

— u;&,(n-Ti£) = — 2^ m-m 0 )exp(27rz — ) (1) 

?n— 0 



32. (New) The method according to claim 19 for accurate pitch estimation, wherein 
the windowed signal is a sound having a pitch and the method further comprises 
accurately estimating said pitch based on the computed frequencies and complex 
amplitudes. 

33. (New) The method according to claim 19, wherein the method is applied to 
determine arbitrary sample rate conversion. 

34. (New) The method according to claim 19, wherein the windowed signal is a 
sound and wherein noise residual, the amplitudes and the frequencies are encoded in a 
bitstream which is stored, broadcasted or transmitted at a sender side of a 
parametric/sinusoidal audio coder, and a receiver decodes the bitstream back to the 
parameters and synthesizes the sound 

35. (New) The method according to claim 19 for audio effects whereby noise r m the 
amplitudes A and the frequencies 6) are manipulated by an effects processor yielding 
C A * and ^>*and synthesized with these modified parameters. 

36. (New) The method according to claim 19 for source separation, whereby 
sinusoidal components originating from the same sound source are grouped and 
synthesized separately. 

37. (New) The method according to claim 19 for automated annotation and 
transcription whereby the signal is segmented according to the values of the amplitudes 
and frequencies. 
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